Molecular dynamics simulation of the human estrogen receptor alpha: contribution to the pharmacophore of the agonists

Human estrogen receptor alpha (ERα) is one of the most studied targets for in silico screening of bioactive compounds. The estrogenic activity of a vast number of chemicals has been studied for their potentially adverse effects on the hormone regulation of the endocrine system. The commonly accepted presentation of the ERα agonist pharmacophore includes terminal phenolic groups and a hydrophobic core with a rigid framework [1]. In this study we report on molecular dynamic (MD) simulations of ERα to get a deeper structural insight into the agonist-receptor interactions and the pharmacophore pattern of compounds with agonistic activity. We rely on a crystallographic structure of a complex of ERα with an agonist of picomolar affinity (PDB ID 2P15, http://www.rcsb.org/pdb/explore/ explore.do?structureId=2P15). As the X-ray structure has mutation in the key structural element for ERα agonistic activity (helix 12, Y537S), a series of MD simulations have been performed on the wild type receptor to prove the stability of the agonist-receptor interactions. The results suggest that the pharmacophore of compounds with ERα agonistic activity can be extended by a feature that occupies a free hydrophobic region of the binding pocket. The results imply also that MD simulations are a powerful in silico tool for both protein dynamics and structure investigation, especially when mutations are available that can potentially disturb the protein structure and functions. The work is supported by the National Science Fund of Bulgaria (grant DTK 02-58).